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ABSTRACT 

We present a compilation of measurements of the stellar mass density as a function of 
redshift. Using this stellar mass history we obtain a star formation history and com- 
pare it to the instantaneous star formation history. For z < 0.7 there is good agreement 
between the two star formation histories. At higher redshifts the instantaneous indi- 
cators suggest star formation rates larger than that implied by the evolution of the 
stellar mass density. This discrepancy peaks at z = 3 where instantaneous indicators 
suggest a star formation rate around 0.6 dex higher than those of the best fit to the 
stellar mass history. We discuss a variety of explanations for this inconsistency, such as 
inaccurate dust extinction corrections, incorrect measurements of stellar masses and 
a possible evolution of the stellar initial mass function. 
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■ 1 INTRODUCTION 

Much contemporary research in extragalactic astronomy has 
, revolved around the determination of the instantaneous cos- 
• mic star formation history (SFH, Madau et al. 1996; Lilly at 
al. 1996). However, measuring this quantity from observa- 
, tions requires a number of assumptions, with the form of the 
dust obscuration corrections and stellar initial mass function 
I (IMF, see Kroupa 2007a for a recent overview) being among 
, the most important. 

Integration of the instantaneous star formation history 
over redshift, making appropriate corrections for stellar evo- 
lution processes, yields the current stellar mass density. This 
quantity can be independently measured, typically using ex- 
tensive galaxy surveys such as the 2dFGRS or SDSS, often 
combined with near infrared (NIR) measurements. Numer- 
ous studies have attempted comparisons of these quantities. 
Madau, Pozzetti & Dickinson (1998), Cole et. al (2001), 
Fontana et al. (2004) and Arnouts et al. (2007) all found 
good agreement between the SFH with a low dust content 
and measured values of the stellar mass density. On the other 
hand there have been a number of studies (Eke et al. 2005; 
Hopkins & Beacom 2006, hereafter HB06) which claim that 
the instantaneous SFH overpredicts the low redshift stellar 
mass density. We attempt to investigate this possible dis- 
crepancy using a compilation of the most up to date mea- 
surements of the stellar mass density history (SMH, the fi*- 
redshift relation). This relation is intricately connected to 



the instantaneous star formation history, but in this context 
it has some important advantages. The principal advantage 
is that estimates of stellar mass typically probe a range of 
the stellar mass function that is somewhat more represen- 
tative of the stellar mass whereas instantaneous indicators 
probe only the most massive stars. Furthermore instanta- 
neous measurements can be subject to a greater uncertainty 
introduced by the effects of dust obscuration. 

In §2 we present a compilation of both low and high red- 
shift measurements of the stellar mass density. Using these 
values, in §3 we derive a best fitting star formation history. 
We then compare this estimate in §4 to other estimates of 
the star formation history and highlight any discrepancies. 
Finally, in §5 we present a discussion of our results and end 
in §6 with a summary. Throughout this work, we assume a 
flat A CDM cosmology with S7a = 0.7, ^matter ~ 0.3 and 
Ho = 70 kms"^ Mpc'^. 



2 THE EVOLUTION OF STELLAR MASS 

The stellar mass density at any redshift is obtained by the 
integral of the galaxy stellar mass function ($, GSMF) 

poo 

p,= I M^{M)dM, (1) 
Jo 

with the GSMF defined such that $(M)dM is the density 
of galaxies with masses between M and M -\- dM. 
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Table 1. Summary of recent measurements of the local stellar mass density (top group) and at higher redsliifts (bottom group), rcdshift 
range, observed mass function range, Salpeter converted value of f2* (with original values in parenthesis) and conversion factor to 
our IMF. "Adopts a diet Salpeter IMF. ''Is independent of an IMF. Adopts a Chabricr IMF. ''Adopts an IMF with a mass range 
of 0.1 — 125Mq. Adopts an alternative cosmology - h = 0.71, = 0.27, f2yv = 0.734. ^ Adopts an alternative cosmology - h = 
0.65, Urn = 0.35, = 0.65. ^Published results only include galcixies with the MF range, the authors estimate they obtain 80% of 
the mass however, we roughly correct by increasing the estimates by 20%. ''Original mass estimates (in square brackets) include only 
galaxies with Ly-'* > 1.4 x lO^^h^gLQ. The authors state that according to the SDSS luminosity function paxameteres they lose 46% 
of light at 2 = 0. At z = 2.8 due to brightening the authors estimate this becomes 30%, thus we correct according to these numbers 
interpolating inbetween.' Values in square brackets are for the Zq, one component ^'(t) model fit with 68% random errors. Utilised 
values and uncertainties are based on the average of model fits using different metallicities and star formation history models. ^ Quoted 
uncertainties are estimates as uncertainties in Foiitana et al. (2006) are quoted only for the stellar mass density in the observed range 
(original values and uncertainties shown in square brackets). * Quoted values are the average of the different methods and samples used 
in the study and errors are based on the scatter between the these methods and samples. Statistical errors are generally < 10% for each 
sample and method. 
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2.1 Stellar Mass Density Measurements at z = 

The most prevalent method of estimating stellar masses is to 
obtain the closest fit between observed galaxy spectral en- 
ergy distributions (SEDs) and a library of template SEDs. 
Template SEDs are typically have a range of star formation 
histories, metallicity distributions and dust content. These 
are typically created using a population synthesis model, a 
code which generates the SED of a stellar population by 
combining a stellar evolution prescription with a library of 
stellar spectra. In addition to requiring a star formation 
history, metallicity distribution and dust content popula- 
tion synthesis models also require the assumption of a spe- 
cific form of the initial mass function (IMF). Variations be- 
tween alternative IMFs can for a given star formation his- 
tory produce different SEDs, and thus different recovered 
masses making this technique then IMF dependent. Com- 
paring results from different studies can then only be ac- 
curately achieved by converting to the same IMF. Further- 
more, different population synthesis models often use dif- 
ferent evolution prescriptions and spectral libraries. This in 
turn can produce variations in the SEDs (Bruzual & Chariot 
2003) of template galaxies and thus a contrasting recovered 
mass (Panter et al. 2007; Pozzetti et al. 2007). 

An early major study using this technique is that of 
Cole et al. (2001) who use galaxy spectra from the 2dFGRS 
and the NIR luminosities from two Micron All Sky Survey 
(2MASS). They find Q^z = 0) = 0.0041 ± 0.0006, where 
is the stellar mass density in units of the critical density as- 
suming the Salpeter IMF. Eke et al. (2005) use an updated 
version of these catalogues to obtain f2, = 0.0033 ± 0.00014 
similarly assuming a Salpeter IMF but using a more gen- 
eralised set of template star formation histories. Bell et al. 
(2003) also use the 2MASS catalogue but instead of the 2dF- 
GRS use the early data release of the SDSS. They obtain a 
stellar mass density of = 0.0028 ± 0.00086, using a mod- 
ified Salpeter IMF (diet Salpeter, Bell & de Jong 2001), a 
Salpeter IMF with a lower limit of M = O.I88M0 motivated 
by numerous observations which suggest the Salpeter IMF 
overpredicts the number of low mass stars (Kroupa 2007a) . 
Panter, Heavens & Jimenez (2004) use the SDSS data re- 
lease 1 obtaining a result of n. = 0.0034 ± 0.00011 for the 
Salpeter IMF. Driver et al. (2007) adopt the diet Salpeter 
IMF and using a survey of galaxies in the B-band obtain a 
result of = 0.0054 ± 0.0008. 

Sampson et al. (2007) determine the local stellar mass 
independent of the IMF, updating the work of Read & 
Trentham (2005). This is achieved by constructing the lo- 
cal typo-specific i?-band luminosity function from a vari- 
ety of local galaxy catalogues. Type and luminosity spe- 
cific mass-to-light ratios obtained from dynamical and lens- 
ing measurements, which do not depend strongly on the 
IMF are then used to determine the GSMF. Integration 
over this GSMF yields fi, = 0.0017 ± 0.0008. Although 
this technique is independent of the IMF it does suffer a 
variety of other limiting assumptions, principally the char- 
acteristics of the dark matter component. In this paper 
we choose to adopt an IMF consistent with recent obser- 
vations of the local IMF (Kroupa 2007a). We define the 
IMF a power law with an index of a = 1.0 over the range 
0.1 < M/Mq < 0.5 and the Salpeter index (a = 2.35) over 
the range 0.5 < M/Mq < 100. Converting estimates using 



alternative IMFs to this IMF is achieved using a population 
synthesis code such as the code Bruzual & Chariot (2003) or 
PEGASE (Fioc & Rocca-Volmerange 1997). Because of the 
ease of introducing alternative IMFs that it allows we adopt 
the PEGASE.2 code. We find that the choice of this IMF 
implies significantly smaller stellar masses than the Salpeter 
IMF, that is Mass(Our IMF) ~ Mass(Salpeter) -0.23 dex 
at ^; = 0. This conversion factor mildly evolves increasing to 
Mass(Our IMF)=Mass(Salpeter)-0.21 dex a.t z = 4. 

2.2 Stellar Mass Density at Higher Redshifts 

Similar techniques as those utilised at z ^ have recently 
been applied to surveys of the higher rcdsliift galaxy pop- 
ulation. The main difference is that surveys at higher red- 
shift often have far fewer galaxies than those locally (such 
as SDSS and 2dFGRS) and thus may be incomplete at the 
extremes stellar mass function. To account for this the mass 
function is often extrapolated beyond the observed range to 
cover a range of typically 10^ < M/Mq < 10^^. This is most 
often achieved by fitting a Schechter (1976) function. 

A compilation of both local and high redshift studies 
is presented in table 1. In this table measurements arc con- 
verted to the Salpeter IMF with the additional conversion 
to our IMF explicitly stated. The measurements converted 
to our IMF arc shown in figure 1. The solid lino in figure 1 
was generated by binning each measurement weighted by its 
uncertainty. A parametric form, p*{z) = a x e~^^ is also 
shown (dashed line) with best fitting parameters a = 0.0023, 
h = 0.68 and c = 1.2. Although this parametric fit works 
well up to intermediate redshifts, at a ~ 4 it is significantly 
smaller than the data. 



3 RECONSTRUCTING THE STAR 
FORMATION HISTORY 

The stellar mass density history p*(t) can be expressed as 
the integral of the star formation history p»(t) corrected 
for the effects of mass loss through stellar evolution pro- 
cesses such as supernovae and stellar winds (Woosley & 
Weaver 1995, Renzini & Voli 1981). Often this is expressed 
as p*{t) = (1 — fi) /q Ptdt' , whore R represents the frac- 
tion of material returned to the ISM. In reality this is an 
approximation as i? is a product of the previous star forma- 
tion history. For example a recent burst of star formation 
would have returned only small amount material where as 
a burst at higher redshift would have returned significantly 
more material by the current epoch. To account for this, the 
stellar mass history can be expressed more accurately as 

p.(t)= / p4t'K^-Mt-t'\)dt', (2) 

Jo 

where fr[t — t'] is the fraction of stellar mass created at t' 
that has been returned to the ISM by t. This quantity can be 
calculated by considering the mass evolution of an instanta/- 
neous burst of star formation using either population syn- 
thesis codes or analytical formulae derived to fit the initial- 
mass/final-mass relations from observations such as those of 
Woosley & Weaver (1995) and Renzini & VoU (1981). We 
obtain this quantity using the PEGASE code to model the 
mass loss characteristics after a burst of star formation. We 
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Figure 1. Stellar mass as a function of redshift converted to our IMF and cosmology. The solid blue line is created by a simple weighted 
bining procedure and the dashed grey line is a best fit with the parameterisation p*{z) = a X e"''^" and a = 0.0023, b = 0.68 and 
c = 1.2. The dotted red line is the prediction of HB06 generated by integrating the observed instantaneous star formation history. 



check the results of this by comparing to both the analytic 
formulae of Hurley, Pols & Tout (2000) as well the alterna- 
tive population synthesis code of Bruzual & Chariot (2003) . 
In the former case we find good agreement for /,. for our 
IMF. Similarly we find that both PEGASE and Bruzual & 
Chariot (2003) provide similar results for both the Salpeter 
and Chabrier IMF. 

Inverting equation 2 can be used to determine the star 
formation history from the observed evolution of stellar 
mass. This yields 

p,{t) = dp,{t)/dt^ dp,.r{t)/dt, (3) 

where dp* {t)/dt is the time derivative of the observed stellar 
mass history (i.e. that presented in figure 1) and p*-r{t) is 
the rate at which material is returned to the ISM. This is 
a function of the previous star formation history and the 
returned function introduced above, 

p,At)= f P4t')ifr[t ^ t'])dt' . (4) 

Jo 

Using the compilation presented in §2 and the formal- 
ism presented in equations 3 and 4 a best fit star for- 
mation history along with la and 3a uncertainty regions 
was generated. For future ease of comparison we follow the 
lead of other authors and express our star formation his- 
tory by the parameterisation of Cole et al. (2001), p, — 
[a + bz)h/[l + (z/c)'']. We find best fitting parameters of 
a = 0.014, b = 0.11, c = 1.4, and d = 2.2 for h = 0.7. This 
best fit star formation history, and the associated la and 3a 
uncertainty regions are shown in figure 2. 



4 COMPARISON WITH OTHER INDICATORS 
OF STAR FORMATION 

4.1 Instantaneous Indicators of Star Formation 

Star formation rates are most often calculated from instan- 
taneous indicators (see Kennicutt 1998 or Calzetti 2007 for 
an overview). This is typically emission associated with on- 
going star formation. Since very massive stars have lifetimes 
which are short compared to typical star formation event 
timescales, emission associated with these stars is used as 
an indicator of their formation rate. These young, massive 
stars completely dominate the integrated UV emission of a 
galaxy, which is thus used as an indicator for their pres- 
ence. Emission from these stars, particularly that which is 
shortward of the Lyman limit is reprocessed by hydrogen to 
produce nebular lines such as Ha and H/3. These along with 
other indicators such as forbidden lines and infrared emis- 
sion (in the case of heavily obscured star formation regions) 
can be used as tracers of the formation of massive stars. 
Extrapolating down the IMF then yields the total star for- 
mation rate. Using the compilation of HB06 (an updated 
version of Hopkins 2004) we compare a range of instanta- 
neous measurements to our star formation history, shown 
in figure 2. The measurements are corrected for the effect 
of dust attenuation using a common obscuration correction 
(as in HB06). As detailed in Hopkins (2004), the common 
obscuration correction corrects emission-line measurements 
using Auct = 1.0 and the Cardelli et al. (1989) Galactic ob- 
scuration curve, and corrects continuum UV measurements 
assuming Av^t„ ~ 0.52 and the Calzetti et al. (2000) star- 
burst obscuration curve. We convert from a Salpeter IMF 
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Figure 2. Comparison of the star formation history inferred from the stellar mass density history compared to other measurements 
of star formation rates. The star formation history inferred from the evolution of stellar mass is shown by the lir and Sct uncertainty 
regions (dark and light grey shaded areas respectively). The dark solid line is the parameterised best fit to our star formation history 
discussed in the text. The lower panel displays the ratio of this best fit SFH to the other measurements. Fossil History results are from 
Panter (2006) (dark green). Instantaneous indicators of star formation are open diamonds colour coded by type and are dust corrected 
as laid out by Hopkins & Beacom (2006). [OH] lines (blue): Teplitz et al. (2003); Gallego et al. (2002); Hogg et al. (1998); Hammer et al. 
(1997). Radio/Sub- mm/FIR (red): Flores et al. (1999); Barger, Cowie & Richards (2000); Condon, Cotton, & Broderick (2002); Sadler et 
al. (2002); Serjeant, Gruppioni, & Oliver (2002); Machalski & Godlowski (2000); Haarsma et al. (2000); Condon (1989); Pcroz-Gonzalez 
et al. (2005). UV indicators (purple): Giavalisco et al. (2004); Wilson et al. (2002); Massarotti, lovino, & Buzzoni (2001); Sullivan et 
al. (2000); Steidel et al. (1999); Cowie, Songaila , & Barger (1999); Treyer et al. (1998); Connolly et al. (1997); Lilly et al. (1996); 
Schiminovich et al (2005); Wolf et al (2003). Xray (grey) Georgakakis et al. (2003). SDSS data (yellow) Baldry et al (2005). Ha and H/3 
(green): Pettini et al. (1998); Perez-Gonzalez et al. (2003); Tresse et al. (2002); Moorwood et al. (2000); Hopkins, Connolly, & Szalay 
(2000); Sulhvan et al. (2000); Glazebrook et al. (1999); Yan, Windhorst, & Cohen (1999); Tresse & Maddox (1998); Gallego et al. (1995); 
Hanish et al (2006). Supernovae Rates are filled red circles from Dahlen et al. (2004) and Cappellaro et al. (2005) 



to our adopted IMF by scaling the values by —0.15 dex ob- 
tained from the ratio of UV luminosities for a simple burst 
using the PEGASE code. Below z — 0.7 the instantaneous 
star formation rate measurements show very good agree- 
ment with our star formation history. At higher redshifts 
however, we see an increasing systematic deviation, with 
instantaneous indicators suggesting much larger star forma- 
tion rates. This deviation appears to peak z 3 where the best 
fit to instantaneous measurements is approximately 4 times 
larger than the best fit to our star formation history. 



4.2 Core Collapse Supernovae Rates 

A similar estimate of the star formation history can be ob- 
tained from the density rate of core collapse supernovae 
(CCSN; Including Type II, lb and Ic supernova). These are 
produced by very massive stars 8 — 50 Mq, at the end of their 
lifetime (Fryer 1999, Dahlen et al. 2004). Because these stars 



are shortlived compared to typical lengths of star formation 
episodes they are contemporary and can thus be used as a 
probe of the instantaneous SFH in a similar way to the pho- 
tometric estimates described in §4.1. Assuming that all the 
stars in this range undergo a CCSN, the rate is related to 
the star formation rate through (Dahlen et al. 2004; Fryer 
1999; Madau, della Valle, & Panagia 1998), 



JnT M-^{M)dM 

P*{z) = — : — ^ Pccsn{z) 

*(M)dM 



(5) 



For our choice of IMF this becomes p*{z) = 
1QQ{Mq) pccsn{z). Observations of CCSN rates con- 
verted into star formation rate densities in this way are 
shown in figure 2. These measurements are completely 
consistent with both the Icr uncertainty region of our star 
formation history and that of instantaneous indicators. 
CCSN rate measurements are however presently limited to 
less than z < 1. 
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4.3 Star Formation Rates from the Fossil Record 

An alternative estimate of the star formation history can be 
generated from the fossil record of star formation in nearby 
galaxies. This involves determining the distribution of stellar 
ages for individual galaxies through the use of stellar popu- 
lation synthesis models. Heavens ct al. 2004 carried out such 
an analysis on the SDSS data release 1. A similar analysis 
was carried out by Panter et al. (2007) using an updated 
version of the SDSS catalogue (SDSS data release 3) and 
higher resolution model spectra. The results of this updated 
version of the analysis are shown in figure 2. 



5 RESULTS 

For z < 0.7 the la uncertainty region of our SFH is con- 
sistent with best fit instantaneous star formation history 
obtained by HB06. This is consistent with the findings of 
recent study of Bell et al. (2007) (a study included in our 
compilation) who find that the instantaneous star formation 
history (from Hopkins 2004) correctly matches the assembly 
of stellar mass over this rcdsliift range. 

For z > 0.7 the best fit star formation history of HB06 
is consistently higher than our star formation history. At 
z — 3 the best fit of HB06 SFH implies a star formation 
rate around 4 times (0.6 dex) larger than that inferred from 
the stellar mass density. This large deviation at high redshift 
offers an explanation for why the integrated star formation 
history implies a local stellar mass density in excess of that 
measured. 

This leads to two possible conclusions, both of which 
may be correct to some degree. Either stellar mass estimates 
are incorrect or the star formation history of HB06 is over- 
estimated at high redshifts {z > 0.7). The former possibility 
requires that both high and low rcdsliift measurements of 
the stellar mass density are incorrect suggesting a systematic 
underestimation. Such an underestimation could be caused 
by a number of possibilities, principally including incorrect 
extrapolation outside the measured GSMF or erroneous cal- 
culation of stellar masses. The latter possibility could be due 
to the often limited extent of photometric information (how- 
ever Fontanna et al. 2006 investigated the effect of including 
Spitzer infrared information on stellar mass determination 
for high redshift galaxies and found little change) or some 
other systematic effect such as metallicity-age degeneracies 
or the effect of dust. Other possibilities include the presence 
of stellar mass missed from surveys in extended galactic ha- 
los (discussed by Bernstein, Freedman, & Madore 2002ab) 
or as free stars in clusters. 

If the instantaneous SFH has been overestimated at 
high redshifts, this could also resolve the issue. A common 
explanation leading to this conclusion is that the effect of 
dust attenuation has been overestimated in measurements 
such as those presented in HB06. Arnouts et al. (2007) for 
example found that using a milder dust attenuation cor- 
rection stellar mass measurements and predictions from in- 
stantaneous star formation rates can be reconciled up to at 
least z = 1.75. The stellar mass estimates of Arnouts et 
al. (2007) are, however, significantly above the average of 
the compilation, as can be seen in figure 1. This suggests 
that even this mild form of dust correction will still not be 



consistent with the average of the compilation. The dust 
correction applied over 3 < a < 6 by HB06 varies between 
2.4 to 3.7 (depending on rest-frame UV wavelength). Rec- 
onciling the HB06 SFH with that inferred from the stellar 
mass density would require scaling down by almost the same 
factors over this redshift range, implying that there be lit- 
tle or no dust obscuration at all over this redshift range. 
This is inconsistent with the recent results of Ouchi et al. 
(2004), Ando et al. (2005), Chary et al. (2005) and Reddy 
ct al. (2007), which show that at least some samples of high 
redshift UV/optically selected star forming galEixies contain 
significant quantities of dust. As noted by many authors 
(e.g., Afonso ct al. 2003), UV/optical selection is strongly 
biased against dusty systems, so the presence of such obscu- 
ration even in UV-selected samples suggests that the extent 
of the obscuration corrections utilised by HB06 is unlikely 
to be overestimated. It is certainly not overestimated by an 
amount sufficient to fully resolve the discrepancy between 
the observed SFH and that inferred from the stellar mass 
density. 

An alternative solution is that the larger star formation 
rates could be explained by an evolution of the star for- 
mation rate calibration (i.e. the conversion factor for con- 
verting a UV luminosity density in a star formation rate 
density). Such a process could speculatively be driven by an 
evolution of the IMF. Investigation of an environmental or 
temporal evolution of the IMF has been carried out by a 
rmmbor of authors. Resolved studies have found little vari- 
ation of the IMF (Kroupa 2007b) with either metallicity or 
environment, however these have generally been limited to 
local galaxies and as such may not be fully representative of 
the entire galaxy population. The possibility of alternative 
IMFs in galaxies with apparently very large star formation 
rates is also discussed in Nagashima et al. (2005), Baugh et 
al. (2005), Le Delliou et al. (2006), Lacey et al. (2007) who 
suggest a fiat IMF to reconcile a number of observations 
with the predictions of the GALFORM (Cole et al. 2000) 
semi-analytical model. Fardal et al. (2007) also investigates 
the possibility of a different IMF in starburst galaxies to ex- 
plain a discrepancy between the extragalactic background 
light, the instantaneous star formation history and the K- 
band luminosity density. 

Within the context of the deviation we observe, an IMF 
which produces more emission associated with instantaneous 
indicators (such as the UV or Ha luminosity) per unit mass 
created is required at high redshift. Such an IMF is likely 
to be top-heavy (or high mass biased) since instantaneous 
indicators are typically dominated by very massive stars. 

In order to investigate this we introduce a simple com- 
posite IMF which has an implicit time dependence, £,{rn, t) = 
fi(Tn) -|- 6(t)$HMB(m). Here ^i{m) is our normal IMF and 
^HMB(?n) is a simple high-mass-biased IMF consisting of 
Salpeter slope truncated at SMq. Introducing this evolv- 
ing IMF changes both the star formation history implied by 
instantaneous indicators, the fraction of material recycled as 
a function of age and the observed stellar mass density. 

Using a simple model for the evolution of the IMF of 
b{t) = 1 — e~'^^ with T = 0.6 we find significantly increased 
agreement, shown in figure 3. In the top panel the best fit to 
the HB06 star formation history for a non-evolving (dotted 
line) and evolving (solid line) IMF is shown alongside the 
star formation history implied by the evolving IMF corrected 



8 Stephen M. Wilkins, Neil Trentham, & Andrew M. Hopkins 



T 



-0,5 - 




12 3 4 



Figure 3. Top panel - The best fit instantaneous star formation 
liistory to the compilation of HB06 obtained with a universal 
IMF (dotted line) compared with that obtained with an evolv- 
ing IMF (solid line). The shaded regions represent the la and 
3(T uncertainty regions of the star formation history inferred from 
the evolvution of the average stellar mass density assuming an 
evolving IMF. Bottom panel - the observed stellar mass density 
assuming a universal IMF (open circles) and assuming an evolving 
IMF (filled circles) as well as the predictions from the instanta- 
neous star formation history assuming a universal IMF (dotted 
line) and evolving IMF (solid line). 

evolution of stellar mass. In the bottom panel the evolution 
of the stellar mass density predicted from the best fit star 
formation history of HB06 for non-evolving and evolving 
IMFs is also shown. The open and filled points denote the 
measurements of the stellar mass density for a non-evolving 
and evolving IMF respectively. In both panels it is clear that 
the inclusion of this evolving IMF reduces the discrepancy. 

This is a simple and ad hoc model, and many other 
forms of evolving IMF may also reproduce the relationship 
between the instantaneous SFH and the SFH derived from 
the evolution of stellar mass. We use this model here simply 
to illustrate the effect of the increasing high mass bias of an 
IMF toward high redshift, and to show that such an effect is 
sufficient to explain the discrepancies between the instanta- 
neous SFH and the SFH inferred from the SMH. In reality 
it is likely that ^{rn, t) has a complex evolution and this is 
currently being investigated in some detail in ongoing work. 



6 SUMMARY & DISCUSSION 

In this work a compilation of stellar mass density measure- 
ments over the range < 2 < 4 was compiled and converted 
into a modified form of the Salpeter IMF. A simple analyt- 
ical parameterisation of the form p*{z) = a x e~^^ with 
a = 0.0023, b — 0.68 and c — 1.2 was found to fit the obser- 
vations well. Using this compilation a best fitting star forma- 
tion history was the obtained. This star formation history 
is well described by the Cole et al. (2001) parameterisation 



with a = 0.014, b = 0.11, c = 1.4, and d = 2.2. This star for- 
mation history was compared to other indicators of the star 
formation history including instantaneous measures, the ob- 
served rates of core collapse supernovae and the results of re- 
cent analysis of the fossil record from the SDSS. Below a red- 
shift of 0.7 there is good agreement between these estimates 
and the la uncertainty region of our star formation history. 
At progressively higher redshifts, the stellar mass density 
and instantaneous indicator inferred star formation histo- 
ries become inconsistent. Instantaneous measures at 2 = 3 
imply best fit star formation rates 4 times larger than those 
inferred from the stellar mass density. 

There are a number of possible causes of this tension. 
These include principally, uncertainty in the cff'ccts of dust 
on both stellar mass estimates and high redshift star for- 
mation rate estimates as well as modeling uncertainties and 
issues regarding the completeness of the galaxy stellar mass 
function. 

In addition, there are more speculative solutions such 
as an efi'ective temporal evolution of the initial mass func- 
tion. We have identified a simple, non-unique model for an 
evolving IMF that reconciles both the SFH and stellar mass 
history. Other recent evidence for an evolving IMF has been 
explored by Dave (2008), and Van Dokkum (2008), who pro- 
vide different parameterisations. A more vigorous investiga- 
tion of the implications of an evolving IMF, and suitable 
choices of parameterisation, is currently underway. 

Given the importance of both stellar mass density and 
star formation rate density measurements to our under- 
standing of the galax;y formation process it is crucial that 
this discrepancy be resolved. To achieve this, improvements 
need to be made to measurements, and extensions, such as 
an evolving IMF, to galaxy formation models need to be 
implemented. Specifically measurements of the stellar mass 
density can be improved through the use of larger and deeper 
surveys to minimise incompleteness, a deeper understanding 
of the eff'cct of dust attenuation and improvements in both 
the template fitting procedure and the underlying popula- 
tion synthesis models used to create the templates. Improve- 
ments in the star formation history can also be made by 
refining models of dust attenuation and using larger deeper 
surveys. Further constraints on the star formation history 
can also come from improvements in core collapse super- 
novae rates studies, and more detailed observations of the 
extragalactic background light, diffuse supernovae neutrino 
background and the chemical evolution. 
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